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Abstract 

In the previous papers (notably, Y. Kodama, J. Phys. A 37, 11169-11190 (2004), and G. Biondini 
and S. Chakravarty, J. Math. Phys. 47 033514 (2006)), we found a large variety of line-soliton solutions 
of the Kadomtsev-Petviashvili II (KPII) equation. The line-soliton solutions are solitary waves which 
decay exponentially in (x,3')-plane except along certain rays. In this paper, we show that those solutions 
are classified by asymptotic information of the solution as \y\ ^ °°. Our study then unravels some 
interesting relations between the line-soliton classification scheme and classical results in the theory of 
permutations. 
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1 The KPII equation and its line-soliton solutions 

The Kadomtsev-Petviashvili (KP) equation 
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where u = u{x,y,t) and = ±1, describes the evolution of small-amphtude, quasi two-dimensional solitary 
waves in a weakly dispersive medium |[T2]| . The case = — 1 corresponding to positive dispersion is known 
as the KPI equation, whereas the negative dispersion (a^ = 1) case is referred to as the KPII equation. The 
KP equation arises in many physical applications including water waves and plasmas (see e.g. [11 J for a 
review). It is a completely integrable system with remarkably rich mathematical structure which is well- 
documented in several monographs ||Tl[T0l[T5l[T8l|20l. Particularly, it has been known that the solutions of 
the KP equation can be expressed in terms of the x-function |[T0ll25]| . 

u{x,y,t) =2^logx{x,y,t). (1.2) 
In this paper, we consider a class of solutions whose x-function is given by the Wronskian determinant 



z{x,y,t) =Wr(/i,...,/;v) 
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with /('^ = d'f /dx' , and where the functions {/nj^^j is a set of linearly independent solutions of the linear 
system 

dy dx^ ' dt 3x3 • ^^-^^ 

In particular, we investigate the line-soliton solutions of the KPII equation, which are real, non-singular 
solutions localized along certain directions in the (x,j)-plane, and decay exponentially everywhere else. For 
example, a one-soliton solution is obtained by choosing N =\m Eq. (11.31) above, and i{x,y, t) = f{x,y, t) = 
gSi -\-e^-, where 

Qrn{x.,y,t)= kmX + k^y + k^J + 6^,0 (1-5) 
with Qmfl, km tor m = 1,2 are constants, and k\ < /c2- The above choices yield the traveling-wave solution 

u{x,y,t) = {{k2-kifsech^\{Q2-Q\)=^0i-r + O)t), (1.6) 

where r = {x,y). The wave vector k := {lx,ly) = {ki —k2,k\ — k\) and the frequency co satisfy the dispersion 
relation, 

-4(0/^ + /^ + 3/2=0. (1.7) 

The solitary wave given by Eq. (11.61) is localized in the -plane along the line L : 6i =62 whose normal 
has the slope c = ly/lx = ki +k2. The one-soliton solution is characterized by two physical parameters, 
namely, the soliton amplitude a = k2 — ki and the soliton direction c = ki+ ^2- The soliton direction can 
be also expressed as c = tan a, where a is the angle, measured counterclockwise, between the line L and 
the positive y-axis. Conversely, any given choice of amplitude (a > 0) and direction of the soliton gives the 
phase parameters ki and ^2 uniquely d&k\ = ^{c — a) and k2 = \{c + a). Note that when c = (equivalently, 
ki = —k2), the solution in Eq. ( 11.61 ) becomes ^-independent and reduces to the one-soliton solution of the 
Korteweg-de Vries (KdV) equation. 



General line-soliton solutions. Like KdV, the KPII equation also admits multi-soliton solutions which 
can also be constructed via the Wronskian formulation of Eq. ( |1.3b by choosing M phases {6m}m=i defined 
as in Eq. (11.51) with distinct real phase parameters ^1 < ^2 < • • • < and then defining the functions 

M 

fn{x,y,t) = anme^"' , n=l,2,...,N, (1.8) 

m=l 

which give finite dimensional solutions of Eqs. (11.41) . The constant coefficients define the N x M co- 
efficient matrix A := {anm), all of whose N x N minors must be non-negative to ensure that the x-function 
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Figure 1.1: Line-soliton solutions of the KPII equation illustrating different interaction patterns: (a) a 2-soliton 
solution, (b) a partially resonant (3,3)-soliton, (c) a Miles resonance (Y-junction). Here and in all following fig- 
ures, the horizontal and vertical axes are respectively, x and y, and the graphs show contour lines of the solution 
M(x,y,f) = 23^ logx(jc,y,f ) for fixed f. 

l{x,y,t) has no zeros in the (x,3')-plane for all t, so that the corresponding KPII solution u{x,y,t) resulting 
from Eq. (11.21 ) is non-singular. 

However, the multi-soliton solution space of the KPII equation turns out to be much richer than that of 
the (l-i-l)-dimensional KdV equation due to the dependence of the KPII solutions on the additional spatial 
variable y. Asymptotically as y ^ ±00, there exist certain (non-decaying) directions which are invariant in 
t, and along which the solution has the form of a plane wave similar to the one-soliton solution in Eq. (11.61 ). 
These asymptotic solitary wave structures, referred to as asymptotic line-soli tons in Ref. [2J, have varying 
amplitudes and directions depending on M, N and the the values of the phase parameters k\,.. . ,kM- More 
significantly, the number N- of asymptotic line-solitons as j — > —00 is in general different from the number 
N+ of the asymptotic line-solitons as y ^ 00, with = M — N and N+ = N. Such multi-soliton configura- 
tions derived from Eq. (11.81 ) are called (A^_,A'^+) -soli ton solutions of KPII |2]|4l. For example. Figure fTTTT c) 
exhibits a (2, 1) -soli ton solution, also known as the Miles resonance solution lUTI . At the interaction ver- 
tex or Y-junction, the three interacting line-solitons with wave numbers and frequencies cOq (a = 1,2,3) 
satisfy the fundamental three-wave resonance condition 

ki+k2 = k3, (Oi+(02 = (03 (1.9) 

The (A'^_,A'^_|_)-soliton solutions exhibit a variety of time-dependent spatial interaction patterns including 
the formation of intermediate line-solitons in the {x,y)-p\a.ne H [I4l [161. In contrast to these nontrivial 
interactions exhibited by the KPII solitons, the KdV multi-soliton solutions experience only a phase shift 
after colUsion. 

A'^-soliton solutions. When = N+ = N (i.e., when M = 2N) the corresponding solutions consist of the 
same number of asymptotic line-solitons as y — > ±00. If in addition, the direction and amplitude of each of the 
N line solitons asy —00 are pairwise equal to each of the N line-solitons as j — > 00, then the corresponding 
solutions are simply referred to as the A/^-soliton solutions of the KPII equation. It will be evident from the 
discussions in the following sections that each A/^-soliton solution can be regarded as a configuration of N 
interacting asymptotic line-solitons where the amplitude and direction of the ?i-th line-soliton are given by 

<^n = kj„-k„i Cn=ki^,+kj^, n=l,...,N. (1.10) 

Thus the n-th line-soliton is parametrized by a pair (^,^_ , kj^^ ) of distinct phase parameters with 1 < /„ < jn < 
2N, or equivalently, by the index pair [/„, j,,]. 

The Y-junction solution found by Miles lITTl is, perhaps the earliest evidence of resonant structure 
present in the line-soliton solutions of KPII. Subsequently, this solution was reconstructed using differ- 
ent algebraic methods in several earlier works (see e.g., lfT9l [8l 1211'). and more general types of resonant 
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line-soliton solutions of KPII were reported in some recent works including Refs. |[T6l l4ll22l. The properties 
of the general line-soliton solutions were systematically investigated in Ref. [2] where these solutions were 
characterized by developing an asymptotic analysis of the x-function. The special case of the A/^-soliton 
solutions was extensively studied in Ref. [141. In particular, an explicit characterization of the A/^-soliton 
solutions space in terms of Grassmannian Gr{N,2N) was presented for the first time in llT4l . In this article, 
we extend the work in Ref. ||21 by giving a combinatorial description of the general line-soliton solutions 
of KPII, and show that the class of the general line-soliton solutions can be enumerated employing these 
combinatorial properties. We also make remarks on those solutions in terms of the positive Grassmann cells. 
Furthermore, we show that the solution space is divided into dual sub-classes of {N,M —N)- and {M —N,N)- 
soliton solutions under the action of space-time inversion {x,y,t) {—x, —y, —t), and give a combinatorial 
interpretation of this discrete symmetry of the KPII equation. Next we consider the A/^-soliton solutions and 
describe how to construct such solutions starting only from the physical data set of N amplitudes and N 
directions of the associated line-solitons. We show that there exists a one-to-one correspondence between 
A/^-soliton solution space and the set of all fixed-point free involutions of the permutation group of 2N ele- 
ments. We exhibit how the combinatorial properties of the A/^-soliton solutions provide a further refinement 
of the A/^-soliton solution space, and thus recover the results of Ref. [14]. 

2 The KPII x-function and asymptotic line-solitons 

In this section, we investigate the general properties and asymptotic behavior of the x-function associated 
with the general line-soliton solutions of the KPII equation. As before, we consider the Wronskian form of 
the x-function given by Eq. (11.31 ) where the functions {/„}jy^i are linear combinations of exponentials with 
M distinct phases as in Eq. (11.81) . Furthermore, we can assume without loss of generality, that the phase 
parameters are ordered d&k\ <k2< ■■■ < ku- 

2.1 Properties of the x-function 

The Wronskian in Eq. (11.31) can be expressed as 

x{x,y,t)=A&t{A@K), (2.1) 

where A = (a„m) is the N x M coefficient matrix, = diag(e^' , . . . ,e^"), and the M xN matrix K is given 
by K = {k"^^), m = 1,2, ... ,M,n = l,2,...,N. Expanding the determinant in Eq. (12.11 ) using Binet-Cauchy 
formula yields the following explicit form of the x-function 

x{x,y,t) = A(mi,...,m/v) exp[e(mi,...,mAr)] Y[ (^m.-^mj, (2.2) 

l<mi<-<mN<M l<s<r<N 

where Q{m\ ,mi^) := 6,„j +6m2 + • • • + 9mN> ^("^1, . . . ,mi^) is the A/^ x A/^ minor of A obtained from columns 
1 < mi < • • • < mp^ < M, and the product term is the Van der Monde determinant obtained using the rows 
1 < mi < • • • < niN <M of the matrix K in Eq. (12.11 ). The basic properties of the x-function following from 
Eq. (12.21 ) are listed below. 

Property 2.1 

(i) The x-function is a linear combination of real exponentials, where each exponential term contains 
combinations of N out of M distinct phases given by 9(mi, . . . ,m/v). Any given phase combination 
9(mi, . . . ,myv) actually appears in the x-function if and only if the corresponding minor A(mi, . . . ,mN) 
is nonzero. Thus there are at most (^) terms in the x-function. 

(ii) If M = A'^, the corresponding x-function in Eq. (12.21 ) contains only one exponential term which gener- 
ates the trivial solution u{x,y,t) = of KPII via Eq. (11.21 ). Hence, M > N for nontrivial solutions. 
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(iii) If rank (A) < A^, then all the N x N minors of A vanish identically, leading to the trivial case x = 
0. Moreover, for rank(A) = N, if all minors A{mi, . . . ,111^) > then x{x,y,t) > 0,'\/{x,y,t) G R^. 
Therefore, the resulting solution u{x,y,t) of the KPII equation is non-singular. 

(iv) The transformation A — > CA where C G GL{N,M.) corresponds to an overall rescaling x det(C)x, of 
the x-function in Eq. (I2.1I ). which leaves the solution u{x,y,t) invariant. This GL{N,R) freedom can 
be exploited to choose the coefficient matrix A in Eq. (11.81 ) to be in the reduced row-echelon form 
(RREF). 

(v) The transformation A AD, & ^ where D e GL{M,R) leaves the x-function in Eq. dlT]) 
invariant. In particular, a diagonal matrix D with diagonal elements dm > 0,m = I, . . . ,M, leaves the 
functions {/„}^ in Eq. (11.81 ) invariant by simultaneously rescaling the m"* column of A by d,„, and 
shifting the each phase constant in Eq. (11.51) as 6,„ —>■ 6,„,o — log(<im)- 

(vi) If for any given n, m in Eq. (II. 81 ). we take /„ = e^'" such that the n^^ row of A has only one non-zero 
entry = 1, then the minors A(mi, . . . ,171^) = 0,m ^ {mi,. . . ,111^}. The resulting x-function can 
be expressed as x(x,3',f) = e^'"Xo{x,y,t), where Xo{x,y,t) contains at most M — I phases (all but 6^) 
which appear in phase combinations of N — I distinct phases. It is then evident from Eq. (11.21 ) that 
x{x,y,t) and Xo{x,y,t) generate the same solution of KPII. Such x-functions are reducible in the sense 
that they can be effectively obtained from a Wronskian ofN —I functions with M —I distinct phases. 

For the remainder of this paper, we will consider the coefficient matrix A to be in RREF. Furthermore, 
to avoid trivial and reducible cases, and to ensure that the solution u{x,y,t) of KPII resulting from the x- 
f unction in Eq. (12.21) are non-singular, we will impose the following restrictions on the coefficient matrix A: 

Condition 2.2 

(i) Positivity: Rank(A) = N < M and all nonzero minors of A are positive. 

(ii) Irreducibility: Each column of A contains at least one nonzero element, and each row of A contains 
at least one nonzero element in addition to the pivot (first non-zero) entry. 

Remark 2.3 The matrices satisfying Condition I2.2f i) above, are called totally non-negative (TNN) matri- 
ces. The classification of the (A'^_,A'^^)-soliton solutions is then given by the classification of the N x M 
irreducible TNN matrices A in RREF. From a more geometric perspective, each TNN matrix parametrizes a 
unique cell in the TNN Grassmannian Gr^{N,M) (see e.g. [23 1), and the classification of the soliton solu- 
tions corresponds to a further refinement of the Schubert decomposition of Gr{N,M) into TNN Grassmann 
cells (see llT4l for the case M = 2N). The refinement is given by a classification of the coefficient matrix A, 
and the minors A(mi , . . . ,111^) represent the Pliicker coordinates of Gr{N,M). We will discuss the geometric 
structure of this classification in a future communication ||6J. 

2.2 Dominant phase combinations and asymptotic line-solitons 

The spatial structure of the solution u{x,y,t) is determined from the asymptotic behavior of the x-function in 
the (x,3')-plane and for finite values of t. Notice that the x-function in Eq. (12.21 ) associated with an irreducible 
coefficient matrix A is a sum of real exponentials with positive coefficients. If only one phase combination 
6(nii , . . . ,mi^){x,y,t) in the x-function is dominant in a certain region of the (x,y)-plane at a given time, then 
X ~ exp{Q{mi , . . . ,mj^)). Consequently, the solution u{x,y,t) of KPII generated by the x-function (12.21 ) is 
exponentially small at all points in the interior of any dominant region, and is localized at the boundaries 
where a balance exists between at least two dominant phase combinations in the x-function (12.21 ). Such 
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boundary is identified by the equation Q{mi ,mi^) = Q{m\ ,m'j^), which defines a line segment in the 
(;c,j)-plane for each t. Note that this phenomenon also arises for the one-soliton solution (11.61 ). which is 
localized along the line 9i = 62 corresponding to the boundary of the two regions of the (jc,y)-plane where 
61 and 62 dominate. In the one-soliton case, these two regions are simply half -planes, whereas in the general 
case the dominant regions could be bounded or unbounded (see e.g. Figure 12.11) . A detailed analysis of the 
asymptotic behavior of the x-function was carried out in Ref. 121, the main results are summarized below. 

Proposition 2.4 The asymptotic properties of x-function il.lh for finite values oft, and for generic values 
of phase parameters ki,.. . , Icm, are as follows: 

(i) The dominant phase combinations of the x-function in adjacent regions of the {x,y)-plane as y ^ 
±00, contain N — I common phases and differ by only a single phase. The transition 6(/,m2, . . . ,mfj) 1-^ 
6(7,m2, . . . between any two such dominant phase combinations occurs along the line defined by 
[i,j] : 9,- = Qj, i ^ j, where a single phase S; in one dominant phase combination is replaced by a phase Qj. 

(ii) Along the single phase transition line [i,j], the dominant phase balance yields 

x{x,y,t) ~ ^Cje^^j''"'—"''^ (2.3a) 

asymptotically as y or as y ^ —00, where the coefficients Ci, Cj depend on appropriate Van der Monde 
determinants and non-vanishing minors of the coefficient matrix A. The asymptotic behavior of the solution 
along [i,j] is given by 

u{x,y,t) ~ ^{kj-ki)^sechmQj-Qi + bij), (2.3b) 
which defines an asymptotic line-soliton. 

(Hi) The number of the asymptotic line-solitons is invariant in time, and so are their amplitudes and direc- 
tions. In particular, the soliton direction is given by the normal direction of [/, f\, which is Cj j = ki + kj, and 
the soliton amplitude is given by aij = \ki — kj\. 

In view of Proposition 12.41 it is natural to denote an asymptotic line-soliton by the index pair [i,j] which 
labels the asymptotic direction as well as the phase parameters in Eq. (I2.3bl ). The asymptotic properties of 
the x-function reveal that there exists a certain pairing {ki,kj), I < i < j < M between the phase parameters, 
which in turn defines an asymptotic line-soliton. However, Proposition 12.41 does not specify how to deter- 
mine these pairings in a given solutions, or equivalently, which phase combinations are actually dominant 
in a given x-function as \y\ —>■ 00. It is still necessary to identify the particular set of asymptotic line-solitons 
associated with any given x-function of Eq. (12.21 ). For this purpose, we first need a result derived in Ref ll2l 
((Lemma 3.1) regarding the dominant phases. 

Lemma 2.5 (Dominant phase conditions) Along the line [i, j] : 6,- = Qj with i < j, the phases 61 , . . . , 6m 
satisfy the following relations, where 6 := 6, = 6^, below. 

(i) As y ^00, e,„ < 6, Vm G {/+ 1, . . . , j - 1}, and 6^ > 6, Vm € {1, l,j + 1, ... ,M}; 
as y ^ —00, 0^ > 0, Vm G {/+ 1, . . . ,7 — 1}, and 0„, < 0, Vm G {1, . . . , / — 1,7 + 1, . . . ,M}. 

The proof of this Lemma based on the ordering ki < ■ ■ ■ < kM is an easy exercise for the reader. Lemma 
12.51 provides a simple yet useful way to determine the dominant phase combinations along the line [i,j]. 
However, a given phase combination Q{mi, . . . ,m!^) can only be dominant if it is in fact, present in the 
x-function of Eq. ( 12.21 ). i.e., if the corresponding coefficient minor A{mi, . . . ,mi^) 7^ in the x-function. 
Therefore, in order to obtain a complete characterization of the asymptotic line-solitons it is necessary to 
consider the structure of the N x M coefficient matrix A in addition to Lemma 12.51 Each asymptotic line- 
soliton [i,j] of a (A^-jA/^^) -soli ton solution of KPII are uniquely determined by a pair of columns of the 
coefficient matrix A, as prescribed below. Once again, the details can be found in Ref. |21. 
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Proposition 2.6 The {N-,N+)-soliton solution of KPII generated from the x-function in Eq. (12.21 ) has ex- 
actly = N asymptotic line-solitons as y ^ oo and N- = M — N asymptotic line-solitons as y ^ — oo. 
The necessary and sufficient conditions for an index pair [i,f\ to identify an asymptotic line-soliton are 
determined by the ranks of two sub-matrices of A defined below in terms of their column indices 

X[/j]:=[l,2,...,/-l,i + l,...,M] F[/j] :=[/+!,... 

The rank conditions are then as follows: 

(i) Each asymptotic line-soliton as y ^ oo is labeled by a unique index pair [e„,j„] with e„ < j„ where 
{^n}n=i l^b^l the pivot columns of A. Moreover, //'rank(X[e„7„]) := r,„ then 

rn<N-\ and rank(X[f'„j„]|f'„) = mn\i{X[enin]\in) = rank(X[e„j„]|e„,7„) = r„ + 1 . 

(ii) An asymptotic line-soliton as y ^ — oo is labeled by a unique index pair with i„ < gn where 
{gn}^=i^ label the non-pivot columns of A. Moreover, //'rank(F[/„g„]) := s^, then 

Sn<N-l and rank(F[/„g„]|/„) = rsmk{Y[i„gn]\gn) = rank(F[/„gJ|/„,^„) = 5„ + 1 . 

Above, {Z\m,n) denotes the sub-matrix Z of A augmented by the columns m and n of A. 

Given the x-function data, which consist of M distinct phase parameters k\,.. . ,kM and a matrix A satisfying 
Condition 12.21 Propositions 12.41 and l2.6l provide an explicit way to identify all the asymptotic line-solitons 
of the corresponding solution of the KPII equation. This method is illustrated via the examples below. 

Example 2.7 Figure IZTT a) illustrates a (2, l)-soliton Y-junction solution ifTTl describing the resonant in- 
teraction of two line-solitons mentioned in Section [T] (see Figure ll.lf c)). This solution corresponds to 
N = 1, M = 3, and is generated by the x-function and the coefficient matrix A, 

x{x,y,t)=e^' +e^^ + e^\ A=(l 1 l) . 

In this case we know from Proposition 12.61 that the number of asymptotic line-solitons as 3^ ^ oo and as 
y — > —00 are one and two, respectively. Applying the rank conditions from Proposition I2.6f i) to the pivot 
column e\ = 1, we see that for the soliton [l,7i] as j ^ 00, rank(X[l 71]) = since N =\. Hence, ji = 3, so 
that the line-soliton is [1 , 3] , which (according to Proposition 12.41) corresponds to the dominant balance of the 
phases 61 and 63 in the x-function. As 3^ ^ —00, Proposition I2.6f ii) implies that for the non-pivot columns 
g\ = 2 and g2 = 3, we should have rank(F[/i2]) = rank(F[/23]) = as well. Consequently, /i = 1, = 2, 
and the resulting line-solitons [1,2], [2,3] correspond to the dominant balance of the phase pairs (61,62) 
and (62, 63), respectively. Thus, the (;c,y)-plane is partitioned in three disjoint regions where each of the 
phases 61 , 62 and 63 dominates, and the solution is localized along the phase transition lines which mark the 
asymptotic line-solitons. 

Example 2.8 Consider N = 2 and M = 4 corresponding to a (2, 2)-soliton solution as shown in Figure lZTl b). 
and generated by the x-function in Eq. (11.31 ) with 

f^=e'^-e'\ f,=e'^+e'^+e'\ A = \ \ " 

The pivot columns of A are labeled by the indices {ei , ^2} = { 1 , 2}, and the non-pivot columns by the indices 
{.?ii^2} = {3,4}. According to Proposition 12.61 the number of asymptotic line solitons are = = 2. 
They are identified by the index pairs [1, ji], [2,72] as j ^ 00, for some ji > 1 and 72 > 2; and by the index 
pairs [/i,3] , [/2,4] as j — > —00, for some ii < 3 and i2 < 4. We first determine the asymptotic line-solitons 
as J ^ 00 using the rank conditions prescribed in Proposition I2.6f i). For the first pivot column ei = 1, 
starting from j = 2 and then incrementing the value of j by one, we check the rank of each sub-matrix 
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Figure 2. 1 : Dominant phase combinations in different regions of the (x,y)-plane (labeled by the indices in parentheses, 
i.e. (/. j) = 9(', y)) and the asymptotic line-solitons (labeled by the indices in square braces) for two different line- 
soliton solutions: (a) a (2, l)-soliton with {ki,k2,kj) = (— IjO,^) at ? = 0; (b) a (2,2)-soliton with (^i,...,A:4) = 
(-l,-i,i,2)atf =0. 



^'[ly]. Proceeding in this way, we find that the rank conditions are satisfied when j = 3: ^[13] = ( j^). So 
rank(X[13]) = \=N-\. Moreover, rank(X[13]|l) = rank(X[13]|3) = rank(X[13]| 1,3) = 2. Thus, the first 
asymptotic line-soliton as j ^ oo is identified by the index pair [1,3]. For 62 = 2, proceeding in a similar 
manner as above we find that j = 3 does not satisfy the rank conditions (since X [23] has rank 2) but j = 4 
does. Therefore, the asymptotic line-solitons as 3^ ^ 00 are given by the index pairs [1,3] and [2,4]. 

Then we consider the asymptotics for j —00. Starting with the non-pivot column gi = 3, we apply the 
rank conditions in Proposition I2.6f ii) to the column / = 2. Then, we have F[23] = 0, and rank(y [23]|2) = 
rank(F[23]|3) =rank(F[23]|2,3) = 1. Hence, the pair [2,3] identifies an asymptotic line-soliton a.sy —00. 
For g2 = 4, we consider / = 1,2,3 and find that the rank conditions are satisfied only for / = 1. In this case, 
F[14] = (0 0), sorank(F[14]) = 1 =A^- 1 andrank(F[14]|l) =rank(F[14]|4) = rank(F [24] 1 1 , 4) = 2. Thus, 
the index pair [1,4] identifies the other asymptotic hne-soliton as j — > —00. 

The dominant phase regions in the {x,y)-pla.ne for this (2, 2)-soliton solution can be also identified from 
Proposition 12.41 First note that the dominant phase combination along x —>■ —00 is given by 6(1, 2) = 61 + 62 
for finite y. This follows from Eq. (11.51) and the ordering ki < k2 < < /C4 of the phase parameters. As the 
slope of the normal to the transition line 6, = 6y, given by the direction Cij = kj + kj decreases from negative 
;c-axis to positive x-axis for 3^ > 0, the asymptotic line-solitons are sorted clockwise as [2,4] and [1,3] as 
y -^00. Thus the dominant phase combinations associated with single phase transitions (cf. Proposition 12.4b 
as J ^ 00 are given by 

6(1,2)^6(1,4)^6(3,4). 

When y < 0, the soliton direction parameter c,y increases from negative x-axis to positive x-axis. Con- 
sequently, the asymptotic line solitons as j — > — 00 are sorted counter-clockwise as [2,3] and [1,4], and 
determine the dominant phase combinations in the (x,y)-plane for j — > —00 as follows: 

6(1, 2)^6(1, 3)-^6(3, 4). 

The asymptotic line-solitons and dominant phase combinations are shown in Figure ITT] where the phase pa- 
rameters are chosen such that C14 = 1 > C23 = 0. Note that in addition to the unbounded dominant regions, 
there is also a bounded region in the x3'-plane where 6(2,4) is the dominant phase combination. The bound- 
aries of this region is formed by the line-solitons [1,4] and [2,3] as j — > —00, together with the intermediate 
line-soliton [1,2]. 

The above examples illustrate how to identify the asymptotic line solitons and the dominant phase com- 
binations of a (A'^_,A'^^) -soli ton solution of KPII from the x-function in an algorithmic fashion. First, apply 
the rank conditions given in Proposition 12. 6l to each pivot column and to each non-pivot column of the coeffi- 
cient matrix A to identify the asymptotic line-solitons as j — > ±00. Next, note that for a x-function in Eq. (l2.2l) 
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associated with a coefficient matrix A in RREF, the dominant phase combination as — > — oo is uniquely 
given by 9(ei , . . . ,ei\i). This is due to the fact that the phase parameters are ordered as < ^2 < • • • < ^m> 
and the coefficient A{ei,. . . ,eyv) of the term ^^('■i. - .'S'a') being the minor of pivot columns, is lexicographi- 
cally the first non-vanishing minor of A with A{e\ ,6^) = 1. Since the line-solitons are sorted according 
to their direction parameter c,y, the dominant phase combinations in the (x,y)-plane can be then determined 
from Proposition 12.41 starting from the dominant phase combination Q{e\ ,6^) as — > —00. 



2.3 Index pairing and derangements 

We show here that the set of unique index pairings in Proposition l2.6[ {[^n, J«]}^=i ^{[hu8n]}^=i identify- 
ing the asymptotic line-solitons as \y\ —>■ 00, has a combinatorial interpretation. Let [M] := {1,2, . . . ,M} be 
the integer set and recall that {ei,...,eN}U{gi,..., Sm-n} is a disjoint partition of [M] . Define the pairing 
map n : [M] [M] according to Proposition 12. 6f i) & (ii) as 

n{e„)= j„,n = 1,2,..., N, %{gn) = in , n = \,2, . . . ,M -N , (2.4) 

where e„ and g„ are respectively, the pivot and non-pivot indices of the coefficient matrix A. Then one can 
show that the map % : [M] — > [M] is a bijection, that is, 71 is a permutation of the set [M]. It is sufficient to show 
that the image 7i([M]) is a set of distinct elements. First assume the contrary, i.e., suppose 7i(/) = = m 
for two distinct elements 1,1' ^ [M\. Then consider the dominant phase combinations associated with the 
single phase transitions (see Proposition I2.4f i)) starting with the Q{e\ ,eN) asx—^ —00, proceeding clock- 
wise, and finally back to 9(ei , . . . ,eN) after a complete revolution. There are altogether M such single phase 
transitions: N transitions as y ^ 0° corresponding to the pivot indices {ei, . . . ,6^}, and M — N transitions as 
y — > —00 corresponding to the non-pivot indices {gi ,gM-N}- Without loss of generality, if we assume 
that the transition 1 —nn occurs before /' —>■ m, then there must be an intermediate m —nn' transition in be- 
tween those two single phase transitions. Now if 6^ G 6(^1 , . . . ,6^) then the transition m —>^ m' must occur 
before the I ^ m transition can take place. Consequently, the intermediate m ^ m' transition can not occur 
as each transition occurs only once during a complete revolution in the (;c,y)-plane. On the other hand, if 
6m ^ 6(^1 , . . . , gA?) then the m ^ m' transition must occur after the I' ^ m transition, and again the interme- 
diate m—^m' transition can not take place. Either way, we reach a contradiction implying that 71 : [M] —>■ [M] 
is one-to-one, and therefore a bijection, thus proving our claim. Note also that 7i has no fixed point because 
n{e„) > en,n = I, . . . ,N and 7i(g„) < g„,n = I, . . . ,M — N. A permutation 71 with no fixed point is called a 
derangement, and an element I G [M] is called an excedance of 71 if 7l(/) > I. We can summarize the above 
discussions as follows. 



Proposition 2.9 The pairing map K defined by Eq. (12.41 ) is a derangement of [M] with N excedances which 
are given by the pivot indices {e\ e/v} of the coefficient matrix A in RREF. 

In Example 12.71 the pivot index ei = I, and the non-pivot indices gi = 2, g2 = 3 for the coefficient matrix 
A, with M = 3. The soliton pairings are: [1,3] as y — > 00, and [1,2], [2,3] as y — > —00. The corresponding 
pairing map from Eq. (12.41 ) is given by 

'1 2 3' 
.3 1 2, 



71 : 



in the bi-word notation of permutation, and 71 has only one excedance: {1}. In Example 12.81 M = 4, the 
pivot and non-pivot indices are {ei = 1, ^2 = 2} and {^1 =3, g2 = 4}, respectively. The pairing map 

/I 2 3 4 
^=^3 4 2 1 

corresponds to the line-solitons [1,3], [2,4] as y ^ 00, and [2,3], [1,4] as y ^ —00 with excedance set { 1 , 2}. 
The pairing map 7i can also be represented by an open chord diagram associated with permutations (see 
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Figure 2.2: The open chord diagram for K = (46523817). 



Q). For example, shown in Figure [Z2l is the diagram for n = (46523817) (using one-line notation for 
permutations). 

The four upper chords in Figure [2]2] originate from the 4 excedances {1,2,3,6} with the right arrows indi- 
cating the increasing order of the numbers. The upper chords correspond to the asymptotic line-solitons for 
J — > oo with the pairings [1,4], [2,6], [3,5], [6,8]. Similarly, the lower chords with reversed arrows indicate 
the line-solitons for y — > — oo namely, [1,7], [2,4], [3,5], [7,8]. Thus, this chord diagram represents a (4,4)- 
soliton solution of KPII. Note here that the [3, 5] -soli ton appears both iny ±00 which corresponds to a 
2-cycle (35) in the permutation 7i. 

Proposition 12.91 furnishes a combinatorial characterization of the (A^_,A'^+) -soli ton solutions of KPII as 
explained below. 

Definition 2.10 Let S+ := {[£H,7n]}^=i (^nd S- := {[injgn]}^^-!^ denote the index sets labeling the asymp- 
totic line-solitons as y ^ 00 and as y ^ respectively. Then two {N-,N+)-soliton solutions of KPII are 
defined to be in the same equivalence class if their asymptotic line-solitons are labeled by the identical sets 
S± of index pairs, where \S+\ := A'^+ = A'^ and \S- \ := N- = M — N. 

Each equivalence class of (A^_,A'^+)-soliton solutions of KPII is uniquely determined by a derangement 71 
as defined in Proposition 12.91 Therefore, each equivalence class of soliton solutions is also associated with 
a unique open chord diagram with = N upper chords and A'^^ = M — N lower chords. Furthermore, in 
view of Remark 2.3, each derangement of [M] gives a unique parametrization of a TNN Grassmann cell in 
Gr'^ {N,M), denoted by IV (ti), and whose dimension can be computed from the number of crossings of the 
corresponding chord diagram Il27ll23]| . In particular, the dimension of the cell W{n) associated with a TNN 
matrix satisfying Condition l2.2l is given by IS 



(Note that the above definitions of C±(/) are switched from those given in Ref. Q). For the example in Fig- 
ure 12.21 the chord diagram associated to 7i = (46523817) has C+(7l) = 2,C_(7i) = 2, and the corresponding 
Grassmann cell in Gr+ (4,8) has dimension 8=4 + 2 + 2. The top (maximal) dimensional cell in Gr+ (4, 8) 
of dimension 16 = 4x4 corresponds to a chord diagram with the maximum number of crossings, i.e. 
Tif^p = (56781234) = (15)(26)(37)(48) with C+{Ti,op) = C-{Titop) = 6, and dimlV(7i,op) = 16 = 4 + 6 + 6. 
An equivalence class of (A^_,A'^+) -soliton solutions of KPII with a given pairing map 71 can thus be as- 
sociated with a unique TNN Grassmann cell W{ti), and the number of free variables parametrizing the 
(A^_,A'^+) -soli ton solution space is given by dimlV(7i). We next proceed to calculate the total number of 
(A^_,A'^+) -soli ton equivalence classes for given values of M and A'^, and enumerate these equivalence classes 
according to the dimensions of their solution spaces. 



dim W (71) = + C+ (71) + C_ (71) 



where the number of crossings C± (71) in the chord diagram are defined by Q 




(2.5) 
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Let Sm denote the permutation group of [M], and let ©m C Sm be the set of all derangements of [M]. Then 
the derangements can be enumerated according to the number of excedances e{ii) of 7i by the generating 
polynomial 

M-l 

Dm{p) = I = I D^^MP" , M > 1 , 

where the coefficients D^ m denote the number of derangements of [M] with A'^ excedances. The total 
number of derangements of Sm is then given by = ^m(I)- The explicit formula for the derangement 
polynomial Dm{p) is obtained from the exponential generating function ll24]| with Do{p) := 1, 

D{p,z) = £ Dm{p)^^ = — ^ (2.6) 

The first few polynomials are given by Di{p) = 0, D2{p) = p, Dj,{p) = p + p^, D4{p) = p + lp^ + p^. 
Moreover, the polynomials Dm{p) are symmetric, that is, its coefficients satisfy 

Dn,m = Dm-n,m^ A^=1,2,...,M-1. (2.7) 

This and various other properties of the derangement polynomial Dm{p) can be found in Ref. ll24l . Note 
that Eq. ( 12.71 ) is equivalent to the relation 

p"Dm{p-') = DMip), 

which in turn follows from the symmetry 

D{p-\zp) = D{p,z) 

of the exponential generating function, and can be verified directly from Eq. (12.61 ). 

The above formulas give the number Dn^m of equivalence classes for the (A'^_,A'^+) -soli tons of KPII for 
a given M = N- +N= and A'^ = A'^^. Thus, when M = 3, there are only 2 classes of line-soliton solutions 
namely, the (2, l)-soliton as in Example 12.71 and also (l,2)-solitons which are related to the (2, l)-solitons 
by the inversion symmetry: {x,y,t) {—x, —y, —t). When M = 4, there are one type each of the (3, 1)- 
and (l,3)-soliton solutions related via the inversion symmetry, but also 7 distinct types of (2,2)-soliton 
solutions. It is also possible to obtain a further refinement of the total number of soliton equivalence classes 
by introducing a ^-analog of the derangement number D^m namely, 

N{M-N) 
r=N 

where D,-^n^m is the number of derangements of [M] with N excedances, and r — N crossings as defined in 
Eq. (12.51) . Then, Di^ m{(i) is the generating polynomial for the Grassmann cells in Gr'^{N,M) corresponding 
to the irreducible TNN matrices satisfying Condition 12.21 and D,-,n^m is the number of those TNN cells 
of dimension r. The upper limit N{M — N) in the sum gives the dimension of the top cell in Gr^{N,M). 
Equivalently, Dr^i^M then gives the number of (A'^_,A'^4_) -soli ton equivalence classes with Nj^ =N,N- = 
M — N, and with r free parameters. Moreover, the total number of (Af_,A'^+) -soli ton equivalence classes is 
given by DN,M{q = 1) = Dn,m- 

It is interesting to note that DN,MiQ) is related to a ^-analog of the Eulerian number |[27]| . 



n 

+ -/-1 



where [k]^j := l+q + q^ + ... + q'^ Ms the ^-analog of the number k. The polynomial E^^n (q) was recently 
introduced in Ref. |27l where a rank generating function for the cells in Gr^{N,M) was derived by building 
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on the work of Ref. |[23l . It follows from Refs. |[27H 7l that the coefficient of q'' in En^m{q) is the number of 
permutations of [M] with N weak exedances, and whose chord diagrams have r — N crossings as defined in 
Eq. (12.51) . Note that / G [M] is called a weak exedance of a permutation 7i G 5m if 7^(0 ^ The following 
result gives the relation between the polynomials Dn^m{q) and Ei\i m{q)- 

Proposition 2.11 For fixed N+ = N, N- = M—N, the generating polynomial for the (N-,N^)-soliton equiv- 
alence classes according to the dimension of their solution spaces is given by 

j=Q \J J 

where E^^n {4) ci^^ the Eulerian polynomials defined above. 

Proof. Let E {N,M) denote the set of all permutations of [M] with N weak exedances, and let D{N,M) CDm 
denote the derangements of [M] with N exedances. Then the polynomials Epf^Mil) and Dpf^Mil) are given in 
terms of the number of crossings c(7l) of the permutation n as 

%eE(NM) %eD(NM) 

Note that for each n < N — I, an element of E{N,M) can be obtained by adding n fixed points to the 
corresponding element of the derangement set D{N — n,M — n). Then for N' = N — n, M' = M — n, 

SclM]. neD(N' ,M') n=0 V " / 

|5|=;i 

Inverting the above formula, yields the desired result. □ 

Proposition 12. 1 1 1 provides an explicit formula for enumerating the (A/^+jA/^-) -soli ton equivalence classes 
according to the dimensions of the associated Grassmann cells in Gr^{N,M). For example, when M = 4 
and N = 2, Proposition 12. 1 1 [ yields D2^A{q) = q^ +4^^ +2q^. This implies that there are one cell of dimension 
4 (top cell), four cells of dimension 3, and two cells of dimension 2. The Gr+(2,4) case will be discussed 
in the next section. 

It turns out that the line-soliton solutions of KPII possess several other combinatorial properties which 
play significant roles in their classification scheme. Some of these properties were addressed in Ref. 1 14 1 
(see also |3|). In this article, we will present these combinatorial structures underlying the line-soliton 
solutions from an algebraic perspective in Section |4l 



3 (2,2) -soliton solutions 

In this section we study the line-soliton solutions of KPII which admit a pair of asymptotic line-solitons as 
\y\ — > oo. But in general, the pair of line-solitons as j ^ oo differ from those as 3^ ^ —00 in their amplitudes 
and directions. We call these the (2, 2) -soli ton solutions, which include the 2-soliton solutions as well. 
Based on our asymptotic results of Section |2l these solutions are specified by prescribing the distinct phase 
parameters {^1 , . . . , ^4} and the 2x4 coefficient matrix A. The x-function of any (2, 2)-soliton is given by 

x{x,y,t)= £ {ks-kr)A{r,s)e^'+^^ . (3.1) 

l<r<i<4 

where A{r,s) denotes the 2x2 non-negative minors of the matrix A. 
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3.1 Classification of (2, 2) -soliton solutions 

For a given set of phase parameters {^i, . . . ,^4}, all possible equivalence classes of (2,2)-soliton solutions 
can be completely enumerated by the derangements of the index set [4] with 2 excedances. Recall from 
Section 1213 that there are altogether 7 distinct equivalence classes of (2,2)-soliton solutions given by the 
coefficient D2^a of p^, in the derangement polynomial Di^{p). The 7 equivalence classes can be further enu- 
merated by the ^-derangement number D2.4 {q) =q^ + + 2q^ according to the number of free parameters 
spanning the solution space of each equivalence class. The coefficients 1,4,2 in D2^A{q) correspond to the 
number of (2,2)-soliton equivalence classes spanned respectively, by 4,3 and 2 free parameters. This is 
illustrated in Figure l3?T] below by the open chord diagrams associated with the (2,2)-soliton solutions. 




(2413) (3142) (2143) 



Figure 3.1: The open chord diagrams for the seven equivalence classes of (2,2)-soliton solutions. 

The number sequence below each diagram in Figure 13.11 is the one-line notation for the corresponding 
permutation % which also designates a certain TNN cell W{%) in Gr(2,4). The left diagram with 2 crossings 
corresponds to the top cell of dimension 4; the middle four diagrams with 1 crossing in each correspond 
to the dimension 3 cells; while the last two diagrams have no crossing, and correspond to the cells of 
dimension 2 in Gr+(2,4). We verify these combinatorial results below by directly analyzing the (2,2)- 
soliton coefficient matrix A which also represents the TNN cells in Gr+(2,4). 

The (2,2)-soliton solutions arise from one of the following two irreducible coefficient matrices in RREF: 

. (\ -a -b\ . (\ a -b\ ^..^ 

^' = [0 1 c ^J' ^^=U Old) ^^-2^ 

with arbitrary non-negative parameters a,b,c,d. We classify them according to the number of independent 
positive parameters in the coefficient matrix A 1 or in A2. Note that the matrices in Eqs. (13.21 ) identify certain 
Schubert cells of Gr(2,4). A further refinement of those Schubert cells is given by the TNN Grassmann 
cells classified below. 

4 positive parameters: There is one such case, which corresponds to the matrix Ai in Eq. (13.21 ) with 
Ai (34) = be — ad > 0. Note that all six exponential terms are present in the x-function in Eq. (13.11) . From 
Proposition 12.61 one finds that the pair of asymptotic line-soli tons are the same as Ij^I ^ 00 with index pairs 
([1,3], [2,4]). The associated pairing map given by the permutation is 7i = (3412) = (13)(24). Thus, this 
case corresponds to an equivalence class of 2-soliton solutions where each asymptotic line-soliton as 3^ ^ 00 
has identical amplitude and direction to another asymptotic line-soliton as j — > —00. The 2-soliton solutions 
will be discussed further in the next subsection. 

3 positive parameters: There are two possibilities namely, A2 with non-zero parameters a, b, d; and Ai 
with 3 free parameters. Consider Ai first. In tliis case neither b nor c can vanish in Ai, else Ai(34) < 0. 
Hence, a 3-parameter family of the matrix Ai with non-negative maximal minors arises in three possible 
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Figure 3.2: Two different (2,2)-soliton solutions of KPll with phase parameters {ki,. . . ,^4) = (— 1, — 5,2). The 
first solution with K — (3421) at f = — 16 and f = 16 is shown respectively, in (a) & (b); while (c) & (d) represent 
the second solution with K = (4312) at f = — 16 and t = 16. Note the space-time reversal symmetry relating the two 
solutions in (a), (d) and in (b), (c). 

ways namely, a = 0,or d = 0, or Ai (34) = but neither of a,b,c,d is zero. That is, 

/•^ . A -b\ , f\ -a -b\ /I -a -b\ ., , ^ 

W^> = (o 1 c ^J' (")^^ = (o 1 c oj' (-)^' = (o 1 c JwtthA,(34)=0. 

(3.3) 

Case (/) in Eq. (13.31 ) gives rise to a (2,2)-soliton solution with asymptotic line-soliton pairs ([2,4], [1,3]) 
as 3; — > 00, and ([2,3], [1,4]) as y ^ — 00. The associated pairing map is given by 7ii = (3421). The second 
matrix in Eq. (13.31 ) corresponds to another (2,2)-soliton solution whose asymptotic line-soliton pairs are 
([1,4], [2,3]) as J ^ 00, and ([1,3], [2,4]) as j — > —00. These two solutions are related via space-time reversal 
{x,y,t) —>■ {—X, —y, —t) as shown in Figure [3^ Furthermore, the pairing map for the second solution %2 = 
(4312) satisfies 7I2 = Ttf The last case in Eq. (13.31 ) yields a (2, 2)-soliton with asymptotic line soliton pairs 
([2,4], [1,2]) as J — > 00, and ([1,3], [3,4]) as j — > —00, which determine the pairing map 713 = (2413). Like the 
previous two cases, case (///) is also related via a space-time reversal to another (2,2)-soliton solution which 
is associated with the coefficient matrix A2 with 3 positive parameters as in Eq. (13.21 ). Thus, A2 corresponds 
to the asymptotic line-soliton pairs ([1,3], [3,4]) as j ^ 00 and ([2,4], [1,2]) as j — > —00, with the pairing 
map7l4 = (3142) =713"^ 

2 positive parameters: It should be clear from Eq. (13.21 ) that if either Z> = or c = in Ai, it can not satisfy 
Condition 12.21 So, A\ will contain precisely 2 free parameters if and only a = d = 0. Similarly, A2 in 
Eq. (13.21 ) will have 2 free parameters if and only lib = 0. Each of the nonzero parameters b,c in A\, and 
a,d in A2 can be rescaled to unity without loss of generality, by utilizing the gauge freedom described in 
Property 12. If v) in Section [2] The matrices resulting in this way from A\ and A2 are denoted by Ap and Aq 
respectively, and are presented below in Eq. (13.41) . Each coefficient matrix produces an equivalence class of 
2-soliton solutions with asymptotic soliton pairs ([1,4], [2,3]) for Ap, and ([1,2], [3,4]) for Aq, as \y\ —>■ 00. 
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Figure 3.3: Three different two-soliton solutions of KPll with the same phase parameters . . . . , k\) = (—2, — 5,0,1), 
illustrating the three equivalence classes: (a) O-type 2-soliton solution, yielding {01,02) = (—5,1) and (01,02) — 
(|,1); (b) T-type 2-soliton solution, yielding (01,02) = (— 2, j) and (01,02) = (2, j); (c) P-type 2-soUton solution, 
yielding (01,02) ~ (— 1, — 5) and (01,02) — (3, 5). (Note that for this P-type solution C[i 4] < cp.3].) 



3.2 Equivalence classes of 2-soliton solutions 

As noted in the previous subsection, there are three types of 2-soliton solutions. They will be referred to as 
0-, T- and P- types (following the terminology introduced in Ref. 1141 ') below. These are identified by the 
canonical coefficient matrices 

/^l 1 OA . /I -1 -1\ , /I -1\ 

with xi > X2 > in Aj. Up to rescaling of columns, the above matrices were also obtained as special 
subcases of the classification presented in Sectional. Aq describes O-type 2-solitons, with asymptotic 
line-solitons [1,2] and [3,4]; Aj describes T-type resonant 2-solitons with asymptotic line-solitons [1,3] and 
[2,4]; and Ap describes P-type 2-solitons with asymptotic line-solitons [1,4] and [2,3]. A distinctive feature 
of the pairing map 7i of a 2-soliton solution is the fact that % is an involution satisfying 7i = ' . As a 
result, 71 can be expressed as a product of disjoint 2-cycles (see e.g., [51). For the permutation group 54, 
there are only 3 such involutions corresponding to the total number of disjoint partitions of [4] into 2 pairs. 
In cycle notation, these involutions are given by tiq = (12)(34), tit = (13)(24) and Tip = (14)(23) for the 
0-, T- and P-type 2-soliton solutions, respectively. The 2-cycles are also evident from the chord diagrams 
illustrated by Figure [3TT] of the previous section where the extreme left (3412) diagram corresponds to the T- 
type 2-soliton, while the extreme right diagrams (4321) and (2143) respectively, represent the P- and O-type 
2-soliton equivalence classes. 

Figure [331 shows a representative solution for each of the three equivalence classes with the same phase 
parameters k\,...,k\. Note that the O- and P-type solitons interact via an X-junction (ignoring the phase 
shifts), while the T-type solitons interact via four Y-junctions connecting the four asymptotic line-solitons 
to four intermediate segments. Each of these intermediate segments satisfy the nonlinear dispersion rela- 
tion (11.71 ) and at each Y-junction the resonance condition (11.91) is satisfied. Thus each intermediate segment 
is also a line-soliton. For example, in Figure [33l b) the asymptotic line soliton [1,3] (as y —00) forms 
the intermediate line-solitons [1,2] and [2,3] at the bottom left Y-junction. The line-soliton [2,3] connects 
with the asymptotic line-soliton [2,4] (as 3; ^ 00) and the line-soliton [1,2] connects with the asymptotic line 
soliton [2,4] (as y — > —00). Similarly, the asymptotic line soliton [1,3] (as 3; — > 00) forms the intermediate 
line-solitons [1,4] and [3,4] at the top right Y-junction. The Une-soliton [3,4] connects with the asymptotic 
line-soliton [2,4] (as 3; ^ 00) and the line-soliton [1,4] connects with the asymptotic line soliton [2,4] (as 

y ^ _oo). 

An important distinction among the three types of 2-soliton solutions is that they cover different regions 
of the soliton parameter space. Suppose (a\,c\) and (02,02) are the soliton parameters of the asymptotic 
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line-solitons of each type with the same set of distinct phase parameters. Since the phase parameters are 
ordered: < • • • < ^4, the soliton parameters satisfy the following relations which can be easily verified 
using Eqs. (II.IOI) . 

(i) For an 0-type 2-soliton solution C2 > c\ and C2 — c\ > ai + a2- 

(ii) For a T-type 2-soliton solution C2 > ci, and |ai — fl;2| <C2 — c\ <ai+a2. 

(iii) For a P-type 2-soliton solution a2> ai and |c2 — ci | <a2 — ai. 

(iv) (c2 -ci)o > (c2 -ci)t > |c2 -ci|p , («!+ 02)0 < («i + «2)t = («i + «2)p , and 
1^22 — ^^ilo = ^2 — It < {a2 — ai)p. 

Note that for O- and T-type solutions the soliton directions are ordered while for P-type solutions the am- 
plitudes are ordered. Any choice of the soliton parameters {a, ,c,- |a, > 0}?^^ would lead to one of the three 
types of 2-soliton solutions provided that {ci itai, C2 ±a2} are distinct real numbers. Thus, the three types 
of 2-soliton solutions divide the soliton parameter space into disjoint sectors bounded by the hyperplanes 
I C2 — c 1 1 = a 1 + £22 and | C2 — c 1 1 = | a i — a2 1 ■ At each boundary between two disjoint regions of the soliton 
parameter space two of the phase parameters coincide. In such situation, it can be shown (by taking suitable 
limits) that the 2-soliton solution degenerates into a Y-junction li4nl4J. 

Yet another difference is in the phase shifts experienced by the asymptotic line-solitons of each type. 
The position of an asymptotic line-soliton [/, j] is determined by the dominant phase combinations across the 
soliton and from the asymptotic formulas Eqs. (I2.3al ) and (I2.3bl ). The pairs of dominant phase combinations 
across the soliton [i,j] as y ^ oo and as y — > — oo are distinct from each other (see e.g.. Figure [331) . This 
fact gives rise to the phase (position) shift A;y = 5^ — 5,^- for the asymptotic line-soliton [/, j] . Since the 
asymptotic positions 6^ have the same linear dependence on time t, the phase shift A,y is independent of t. 
The dominant phase combinations across an asymptotic line-soliton in each of the three types of 2-soliton 
solutions can be determined from the asymptotic analysis of section[2]and are shown in Figure [331 Then the 
soliton phase shifts computed using the parameters (ai,ci), {a2,C2) and the canonical coefficient matrices 
in Eq. (13.41) are given by 



The phase shifts experienced by the pair of asymptotic line-solitons are of opposite signs in each of the three 
cases. In above, we assume that ci 7^ C2 (i.e., the solitons are not parallel), then A denotes the phase shift of 
the soliton with direction parameter C2 where C2 > ci. For non-resonant (O- and P-type) 2-soliton solutions, 
the phase shift has only one and the same term which depends symmetrically on the soliton parameters. 
However it is easy to verify from the various inequalities mentioned above, among the soliton parameters 
that Aq is always positive while Ap is always negative. For resonant 2-soliton solutions there is an additional 
term that depends on the free parameters xi , X2 of the coefficient matrix Aj in Eq. (13.41 ). In this case, the 
phase shift Aj is not sign-definite unlike the other two cases. 

4 Duality and /V^-soliton solutions 

In the previous section, we described two types of equivalence classes of line-soliton solutions of KPII: 
(i) the (2,2)-solitons where the sets S± (cf. Definition 12.101 ) of asymptotic line-solitons as y ±00 are 
distinct, and (ii) the 2-soliton solutions with 5_ =5+. Furthermore, we noted that there are pairs of distinct 
equivalence classes of (2,2)-soliton solutions related via space-time reversal (cf. Figure [331 ). The general 
(A^^-jA^^-i-) -soliton solutions can be also categorized in a similar fashion according to whether the sets of 
asymptotic line-solitons 5+ and S- are distinct, or if 5_ =5+. In the first case, pairs of equivalence classes 
of solutions are related by space-time reversal, while the latter corresponds to the special case of A/^-soliton 
solutions to be discussed in this section. 



Ao = log 



(ci -C2f - {ai-a2Y 
{c\-C2Y-{ai+a2Y 




(3.5) 
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4.1 Duality of line-solitons 



The KPII equation (11.11) is invariant under the inversion {x,y,t) —>■ {—x,—y,—t). Consequently, if u{x,y,t) 
is a (A'^_,A'^^)-soliton solution of KPII with given sets S± of asymptotic line-solitons and with A'^_ = M — N 
and A'^^ = A'^, then u{—x, —y, —t) is a {N,M — A'^)-soliton solution with reversed sets of asymptotic line- 
solitons. We refer to the solutions u{x,y,t) and u{—x, —y, —t), as well as their respective equivalence classes 
as dual to each other. Let XjMMixjy,!) denote the x-function in Eq. (12.21 ). generating the solution u{x,y,t) 
via Eq. (II. 21 ). then the solution u{—x,—y,—t) will be generated by XNM{—x,—y,—t) as Eq. (11.21 ) remains 
invariant under {x,y,t) {—x, —y, —t). Note however that Xnm{—x, —y, —t) does not have the form given 
by Eq. (12.21 ). but it is indeed possible to construct a certain x-function iM-N.M{x,y,t) from Xn^m{—x, —y, —t) 
that is dual to ijyi^M{x,y,t). We describe below how to construct the function Xm-nm from Xa?m- 

First we obtain a coefficient matrix for the x-function Xm-nm from the N x M matrix A associated with 
Xi\iM(x,y,t). Since A is of full rank, its rows form a basis for a A/^-dimensional subspace W of R.^. Let 
W-^ be the orthogonal complement of W with respect to the standard inner product on R*', and let A' be 
a. {M — N) xM matrix whose rows form a basis for W^. Clearly, A' is not unique, but a particular choice 
for A' is as follows: Suppose the pivot and non-pivot columns of A in RREF are represented by the identity 
matrix 7^ and the A'^ x (M — N) matrix G respectively, then 

A = [In,G]P A' = [-G^ ,Im-n]P , (4.1) 

where G^ is the matrix transpose of G, and f is a M x M permutation matrix satisfying = It can be 
directly verified from Eq. (14.11 ) that AA'^ = 0, which constitute the orthogonality relations among the row 
vectors of A and A'. Moreover, it was shown in Ref. [2J that the A/^ x minors of A and the {M—N)x (M—N) 
minors of the matrix 

B':=A'E, £ :=diag(-l,l,-l,...,±l) 

are related as 

A(mi, . . . = (-1)^ det(P) B\h , ■ ■ ■ Jm-n) , 

where a = M(M+ l)/2 + A^(A'^+ l)/2, and where {mi,m2, ■ ■ ■ ,niN} is the complement of {/i,/2, • • • Jm-n} 
in [M]. The factor (—1)'' det(P) = ±1 which depends only on A, can be rescaled by an orthogonal transfor- 
mation B' — > OB' , det(O) = ±1, so that the maximal minors of the rescaled matrix 

B := OB' = 0[-G^ ,Im-n]PE , (4.2a) 

satisfies the precise complementarity conditions 

A(mi,...,m/v) =B(Zi,...,ZM-iv). (4.2b) 

Hence, if all maximal minors of A are non-negative, then the same holds for B. 

The (M — N) xM matrix B plays the role of a coefficient matrix for the x-function 'ZM-N,M{x,y,t) which 
is related to the function Xa?j^(— a;, —y, —t). Indeed, if we initially set 6„, o = 0, Vm = 1, . . . ,M in Eq. (12.21 ) 
using Property |2. If v). then under the transformation {x,y,t) — > {—x, —y, —t), Eq. (12.21 ) yields, 

Xiv,M(-^, -y, -t) = exp[-e(l, . . . ,M)] x'{x,y,t) . 

Using Eq. (I4.2bl ) and taking the sum over the complementary indices Zi , . . . , Im-n (instead of mi , ... , m/v), 
the function x'{x,y,t) can be expressed as 

x'{x,y,t)= ^ V(mi,...,mA,)B(/i,...,ZM_A?) exp[e(Zi,...,/M-iv)], (4.3) 

l<h<-<lM-N<M 

where 

y(mi,...,m/v) = Y\ {kmr-Ki,) 

l<s<r<N 
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are the Van der Monde coefficients as in Eq. (I2.2I ). Next, if we introduce new phase constants by 



52 ln\kT-k, 



m = \,...,M 



which satisfy the identity 



and make the replacement: 6, 



exp[e/,,o + - 



'm 




x'{x,y,t) 



£ V(Zi,...,/M-iv)fi(/i,---,/M-/v) exp[e(/i,...,ZM-iv)]- (4.4) 



M) 



l<h<-<lM-N<M 



It is clear from Eqs. (11.21 ) and (14.41 ) that the functions Xn,m{—x, —y, —t), x'{x,y,t), and XM-N.M{x,y,t) give 
rise to the same solution u{—x, —y, —t) of KPII. Thus we have the following. 

Proposition 4.1 (i) The equivalence classes of solutions generated by the N xM coefficient matrix A and 
the (M — N) X M matrix B defined by Eq. (I4.2al ). are dual to each other If a (M — N ,N)-soliton solution 
u{x,y,t) of KPII belongs to a certain equivalence class, then its dual equivalence class contains the {N,M — 
N)-soliton solution u{—x, —y, —t). 

( ii) Let {m\ , . . . , mj^} and {Zi , . . . , Im-n} be a disjoint partition of the integer set [M\, then 6(mi , . . . , mj^) is 
a phase combination present in the x-f unction Xn,m if and only ifQ{li,... ,Im-n) is a phase combination in 
the dual x-function Xm-nm- 

( Hi) Ifn G Sm is the pairing map for a given equivalence class, then the pairing map for the dual equivalence 
class is given byn^^. 

Proposition 14. 1 [ establishes a one-to-one correspondence between an equivalence class and its dual. Indeed 
from Eq. (12.71 ) in Section |2]3, we have the following result for the dual equivalent classes. 

Proposition 4.2 Forgiven positive integers M,N, with N <M, the number of of distinct equivalence classes 
of the (M —N,N)-soliton solutions are exactly the same as the number of dual equivalence classes of{N,M — 
N)-soliton solutions. 

Remark 4.3 The open chord diagrams of the pairing maps 71 and ' which correspond to a line-soliton 
equivalence class and its dual aie related to each other via a reflection about the horizontal line together 
with reversing the direction of the chords. This is due to the fact that the exedance set of ' is given by 
{'^{gn)}^Ji' , while the anti-exedance set is given by {i^{en)}^^i - This transformation on the chord diagrams 
can be regarded as the combinatorial analogue of the inversion symmetry {x,y,t) {—x, —y, —t). acting on 
the solutions u{x,y,t) of the KPII line solitons. 

When M = 2N, it follows from Proposition 12.61 that N- = N+ = N, which leads to the (A/^,A/^)-soUton so- 
lutions. In this case the number of asymptotic line-solitons as 3^ ^ oo and as 3^ ^ —00 are the same, but in 
general, the amplitudes and directions of the line-solitons will be different as seen for the (2,2)-soliton ex- 
amples in Figure [l!2l of Section[3l A particularly interesting subclass of the (A'^,A'^) -solitons are the A/^-soliton 
solutions which were introduced in Section [H and which are characterized by identical sets of asymptotic 
line-solitons as \y\ —>■ 00, i.e., S- = 5+. We discuss them next. 

4.2 The A^-soliton solutions 

This special family of line-soliton solutions of KPII consists of equivalence classes that are invariant under 
the inversion symmetry {x,y,t) (— x, —y, —t). That is, each equivalence class is its own dual so that both 
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u{x,y,t) and u{—x,—y,—t) belong to the same equivalence class of solutions. This implies that for each 
A^-soliton solution u(x,y,t), the asymptotic line solitons arise in pairs, where each pair consists of a line 
soliton as J — > oo and another soliton as y ^ — oo, moreover, both line-solitons have identical amplitude and 
direction. Thus, a (A^-jA/^^) -soli ton solution of KPII is a A/^-soliton solution if and only if it is self-dual, i.e., 
if and only if the index sets labeling the asymptotic line-solitons satisfy 5_ =5+. The main features of the 
A/^-soliton solutions follow from from Propositions I2.6IIZ91 14. H and Definition l2. 101 These are listed below. 

Property 4.4 

(i) The x-function of an A'^-soliton solution, denoted by Xn := x^^n, is expressed in terms of 2N distinct 
phase parameters and an N x 2N coefficient matrix A which satisfies Condition 12.21 Then it follows 
from Eq. (I4.2bl ) that the N x N minors of A satisfy the duality conditions: 

A{mu...,mN)=0 ^ A(/i,...,/w) =0, (4.5) 

where the indices {mi, . . . ,mN} and form a disjoint partition of integers {1,2, . . . ,2A^}. 

That is, the phase combination Q{m\ m/v) is present in % if and only if 6(Zi , . . . , In) is. 

(ii) Each A^-soliton solution exactly asymptotic line-solitons as j — > ±00 identified by the same index 
pairs [e„ , g„] with < gn, n = \, . . . ,N. The sets {ei,...,e]\j} and {gi,..., g^} label respectively, the 
pivot and non-pivot columns of the coefficient matrix A. Hence, they form a disjoint partition of the 
integer set [2N] . 

(iii) The amplitude and direction of the asymptotic line soliton [en,gn] are the same as 3^ ^ ±00, and 
are given in terms of the phase parameters as a„ = kg^^ — and c„ = kg^^ + . 

(iv) The pairing maps associated with A/^-soliton solutions are involutions of Sin with no fixed points, 
defined by the set I2N = {ti G S2n\'^^^ = ^7 ^(0 7^ i, V/ G [2A'^]}. Such permutations can be expressed 
as products of disjoint 2-cycles, and theirs chord diagrams are self-dual, i.e. symmetric about 
the horizontal axis (see e.g., the 2-soliton chord diagrams in Section [3] The total number of such 
involutions is given by \I2n\ = {2N - l)!! = l- 3...- (2A^ - 1) t5j|. Hence, there are (2A^ - 1) ! ! distinct 
equivalence classes of A'^-soliton solutions. 

Examples of the A'^-soliton solutions with special choices of the functions {/„}^^i in Eq. (11.81 ) and the 
coefficient matrix A are given below. 

Example 4.5 O-type A'^-soliton solutions. These are the well-known ||9j|T5l multi-soli ton solutions of KPII 
constructed by choosing {fn}^^[ according to 

Mx,y,t)=e^^-"-^+e'''", n = \,...,N. 

The corresponding coefficient matrix is given by 

/I 1 ••• 0\ 
1 1 ••• 



\0 ••• 11/ 



with A'^ pairs of identical columns at positions {In — l,2n}, n = 1, . . . ,N. Thus, there are 2^ nonzero max- 
imal minors of Aq, given by Ao{mi ,171^) = 1, where m„ = 2?i — 1 or m„ = 2n for « = 1, . . . ,N. The A'^ 
asymptotic line-solitons are identified by the index pairs {\2n— 1 , 2?i]}^^j , and the corresponding permuta- 
tion is 71 = (2, 1,4,3, . . . ,2N,2N-l) (or 7i = (12)(34) • • • (2A^,2A^- 1) in the cycle notation). The amplitude 
and direction of the w-th soliton are a„ = — '^2n-i and c,j = ^2n-i +^2n> respectively. Note that the soliton 
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Figure 4.1: Three different 3-soliton solutions of KPII with the same phase parameters {ki,...,k(^) = 
(— 3,-2,0, 1, — |, 2), illustrating the three equivalence classes: (a) O-type, (b) T-type, and (c) P-type 3-soliton so- 
lutions. 

directions are ordered as ci < C2 < ■ ■ ■ < cn due to the ordering of the phase parameters k„. In fact, the soliton 
parameters satisfy the inequalities: c„+i — c„ > a„ + a„+i, n = 1,2,. . . ,N — I. Therefore the asymptotic 
line-solitons for the O-type solutions can not take arbitrary values of amplitude and direction, and thus do 
not cover the entire soliton parameter space as was already noted for the 2-soliton case in Section[3l2. Apart 
from the position shift of each soliton, the interaction gives rise to a pattern of N intersecting lines in the 
(x,3')-plane, as shown in Figure |4~]T a). 

T-type A^-soliton solutions. These solutions are obtained by choosing the functions in Eq. (11.81 ) as 



which yields the coefficient matrix 



2N 



n = \,...,N with f{x,y,t)='£e^ 

m=l 



/ 1 



1 



jN-l 
K2 



1 \ 

kiN 



'^2N / 



In this case all of the NxN minors of the coefficient matrix A are positive, each being equal to a Van der Monde 
determinant. These solutions were investigated in Ref. [4J where it was shown that they also satisfy the 
finite Toda lattice hierarchy. The line-soliton is labeled by the index pair [n,n + N]; that is, its ampli- 
tude and direction are determined by the phase parameters {k„,kn+N)- The corresponding permutation is 
71 = (A^+ 1,A^ + 2,...,2A^, 1,2,...,A^) (or 71 = (1,A^+ I)(2,A^ + 2) • • • (A^,2A^)) whose open chord diagram 
has the maximum number of crossings N{N — l)/2. This implies that the T-type A'^-soliton solution be- 
longs to the top cell of Gr'^{N ,2N). Like the O-type A/^-soliton solutions, these T-type A/^-soliton solutions 
do not cover the whole soliton parameter space. In this case the soliton parameters satisfy the constraints: 
kn-hi —^n\< Cn+\ — Cn < ^n+i +<^n, u = 1 ,2, . . . ,N — I. It was also shown in Ref. |4] that these soliton 
solutions display phenomena of soliton resonance and web structure as shown in Figure l4~]T b). Moreover, 
the intermediate interaction segments are also line-solitons because they satisfy the dispersion relation (11.71 ). 
All of the asymptotic and intermediate line-solitons interact via three-wave resonances. That is, at each in- 
teraction vertex or Y-junction (cf. Figure [TTTJc), the three interacting line-solitons satisfy Miles' resonance 
condition given by Eq. ( 11.91 ). 

P-type A^-soliton solutions. Yet another type of A^-soliton solutions is obtained by prescribing 



fn{x,y,t) = + (-If , n = \,2,...N, 
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in Eq. (11.81) . The coefficient matrix is given by 



/I 
1 



Ap 








*\ 
* 



••• 1 -1 ••• 
\0 ■■■ 1 1 ••• 0/ 



where the asterisk in the {2N — n + 1)* column is equal to (—1)'^^". Like Aq, the matrix Ap also has 
pairs of parallel columns labeled by {(«, 2N — n+ l),n = 1,... ,A^} and 2^ non-vanishing minors, and 
each nonzero minor is 1 . The soli ton is identified by the index pair [n, 2N — n + l], and the corresponding 
permutation is 7i = (2A^,2A^- 1,...,2,1) (or 71 = (1,2A^)(2,2A^- 1) • • • (A^,A^+ 1)). The soliton direction and 
amplitude are given by c„ = kn + k2N-n+i and a„ = k2N-n+i —Ku respectively. Note that the soliton directions 
are not ordered as in the previous two cases. In fact, taking ci = C2 = ■ ■ ■ = cn = yields the reduction to 
solutions of the KdV equation fTTj. But the soliton amplitudes in this case are ordered as ai > a2 > • • • > <^n- 
Moreover, the soliton parameters satisfy the constraints : | c„+ 1 — c„ | <an — a,i+ 1, ?i = l,2,...,A'^— 1. These 
solutions interact non-resonantly, like the 0-type A/^-solitons, i.e., pairwise with an overall phase shift after 
collision (see Figure ITO c)). However, as in the case of the 2-soliton solutions in section |3]2, the pairwise 
phase shifts for P-type solitons is of opposite sign from that of 0-type soliton solutions. 



4.3 Soliton parameters and pairing map 

The most notable differences between the 0-, T- and P-type A/^-soliton solutions in Example |43] are that they 
span different regions of the soliton parameter space and that they exhibit dissimilar interaction patterns 
and phase shifts. However, in addition to these non-resonant and fully resonant solutions, a large family 
of partially resonant solutions exists when N > 2. Thus, the family of A/^-soliton solutions of KPII is much 
larger than previously thought, and their classification is indeed nontrivial. Even forN = 3, Property [4. 4f iv) 
implies that there are 15 inequivalent types of 3-soliton solutions. For increasing values of N, it turns out to 
be a difficult task to classify these solutions according to their coefficient matrices, as was done for = 2 
in Section [3] Instead, a more direct approach is to enumerate the A/^-soliton solutions via the involutions 
hN C Sin, by constructing a representative coefficient matrix A for each A/^-soliton equivalence class starting 
from a pairing map n € hN- In what follows, we describe a slightly modified approach. We start with the set 
of amplitudes and directions of the asymptotic line-solitons as \y\ oo. We first recover the soliton pairings 
from this physical data, then construct the coefficient matrix A from the obtained pairing map. This provides 
a method to algebraically reconstruct the A/^-soliton solution unique up to space-time translations, starting 
simply from the physical soliton parameters {(an,c„)}^^j. 

We begin with the following definition of the A/^-soliton parameter space. 

Definition 4.6 An N-tuple of pairs pj^ := {(fl„,c„)|a„ > 0}^^j C M^^, of amplitudes and directions for the 
asymptotic line-solitons associated with an N -soliton solution, is defined to be admissible if it yields the set 
'■= {{Ki ■,kn)\kn = ^(cn ±a«)}^=i where the 2N phase parameters {kn}n=\ distinct. The set Sol{N) 
of all admissible N -tuples of amplitudes and directions will be referred to as the N -soliton parameter space. 



Note that the pairs {an,Cn),n = l,...,N in the set p^ are unordered. For example, p2 = {(1,2), (^,1)} 
and p2 = {(^,1),(1,2)} represent the same 2-soliton solution. Similarly, IT/v also consists of unordered 
pairs {k^,k^),n = l,...,N. However, since the parameters {kf}^^i are distinct, they can be sorted in 
increasing order into an ordered set K2N = {^1 ,^2, • • • ,^2^}- Hence, YIn forms a partition of K2N into N 
distinct pairs. The positions of each pair {k^ ,k^) e Il^v can be uniquely identified within K2N by an ordered 
pair of indices [/„, j„] such that /„ < 7',,. That is, k^ = ki^^ and k^ = kj^^. It is precisely this identification 
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that induces a correspondence between each Yl^f and a pairing map n e I2N C 52^, the latter representing a 
disjoint partition of [2N] into N distinct pairs (see Property l4.4f iv)). It is then clear from Definition 14 . 61 that 
this correspondence also extends between each pa? G Sol{N) and a 7i G l2N- 

Example 4.7 Consider a 3-soliton parameter set pi, = {{ai,ci) = (1, — 3), (aijQ) = (§, 5)) (^^35^3) = 
(|, |)}, and construct the set 113 = {^{cn ^<^n)}n=i- The set 773 is admissible because the corresponding set 
113 = {(—2,-1), (—5, 1), (2,2)} contains six distinct phase parameters. Sorting these pai^ameters in increas- 
ing order yields K(, = {-2,-1,-^,^,1,2} = {^1,^2,^3,^4,^5,^6}- Then, 113 = {(^1,^2), (^4,^6), (^3,%)} 
which gives the correspondence 773 ~ 113 (12) (35) (46) = 71. 

Note however that the correspondence between the soliton parameter space Sol{N) and I2N is not one-to- 
one because distinct elements of Sol{N) associated with solutions in the same equivalence class give rise to 
identical pairing. In this situation, the corresponding sets Yl^ are distinct but after sorting, their elements 
are ordered in identical fashion into the respective sets K2n- Thus, the soliton parameter space Sol{N) is 
partitioned into disjoint sectors. Each sector corresponds to an equivalence class of solutions, distinguished 
by an element 71 S l2N^ equivalently, by the set {[en,gn]}^^^i labeling the A'^ asymptotic line-solitons. The 
total number of such disjoint sectors of Sol{N) equals the cardinality |/2Ar| = {2N — 1)!!. 

Once a pairing map 7i G I2N is derived from a given A/^-soliton parameter set p^, it is then possible to 
construct a coefficient matrix A satisfying Condition 12.21 Clearly, the pivot columns of A will be labeled by 
the excedance set {ei ,ei\i} of 7i, and the non-pivot columns are labeled by {g\ ,gN}, where 7l(e„) = 
gn,n = I, . . . ,N. The explicit form of A will be determined by using the rank conditions in Proposition 12.61 
Recall that in Examples 12.71 and 12.81 we demonstrated how to apply the results of Proposition 12.61 to a 
given coefficient matrix A, and obtain the set of index pairs identifying the asymptotic line-solitons. In 
the examples below, we will illustrate the reverse construction. In other words, we will show that the rank 
conditions of Proposition I2.6l are also sufficient to construct a coefficient matrix A from a given pairing map 
71 associated with the A'^-soliton solutions. 

Example 4.8 We outline the construction of a coefficient matrix A associated with the 3-soliton pairings 
{[1,2], [3,5], [4,6]} found in Example 14.71 The construction proceeds in several steps. 

Step 1. It follows from Property l4.4f ii) that the pivot and non-pivot columns of the 3x6 matrix in RREF 
are labeled by {^1,^2, ^3} = {1,3,4}, and {^1,^2, ^3} = {2,5,6}, respectively. So, the general form of A 
satisfying Condition I2.2l is 

/I z Vl Wl \ 
A = 1 -V2 -W2 , 
yO 1 V3 W3 / 

where z,{v,-,w,}?^j ai^e non-negative numbers to be determined. Note that each unknown entry can be 
expressed as certain maximal minors of A, (e.g., A(234) = z, etc). The negative signs in the second row are 
included so that the maximal minors of A satisfy the non-negativity Condition 12. 2f i). 

Step 2. In order to obtain further information about A, one needs to apply the rank conditions in Proposi- 
tion l2.6l to the sub-matrices X[ij] and Y[ij] associated with each line soliton 

If we start with the line-soliton [1,2], and consider the sub matrixF[12] =0, we find that rank(F[12]) =0. 
Then according to Proposition 12. 6f ii). rank(F [12] 1 1 . 2) = 1 , so that columns 1 and 2 are proportional. Hence, 
Zy^O. The sub matrix X [12] consists of columns 3, . . . , 6 of A above. From Proposition 12. 6f i). it follows that 
rank(X[12]) < A'^ — 1 = 2, since N = 3. Then any maximal minor of A consisting of the columns of ^[12] 
will vanish. In particular, we find that A(345) = vi = and A(346) = wi = 0, so that 

/I z \ 

A = 1 -V2 -W2 , 
\0 1 V3 W3 / 
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Step 3. Next we consider the index pair [3,5]. As the sub matrix F[35] = (0,0, 1)^ is of rank 1, we obtain 
ranlc(y[35]|5) = 2 from Proposition l2.6f ii). Therefore, columns 4 and 5 of A are linearly independent. Then 
V2 7^ 0. The sub-matrix X[35] consisting of the columns 1,2 and 6 of A, is of rank 2, and its columns space 
is spanned by the linearly independent columns (1,0,0)^ and (0, — W2, wa)^. Then from Proposition I2.6f ii). 
rank(X[35]|3) = 3 ^A(136) = >V3 7^0, andrank(X[35]|5) = 3 ^A(156) =V3W2-V2W3 /O. Moreover, the 
non-negativity of all maximal minors of A requires that z,V2,h'3 and A(156) must be positive. In particular, 
A(156) = V3W2 — V2W3 > implies that V3 / 0, W2 / 0, so these are also positive. Thus, the matrix A is 
parametrized by 5 positive parameters which can be chosen as follows: ti = zji = ^2,^3 =vj,t4 = wsJs = 
A(156) = V3W2 — V2H'3. It can be directly verified that all non-zero maximal minors of A are polynomials in 
ti,...,t5 with positive coefficients. 

Thus, we have constructed a 5-parameter family of coefficient matrix A that will generate the asymptotic 
line-solitons [1,2], [3,5], [4,6] for any choice of the positive parameter values t\,...,t5, via Proposition 12.61 
Note that the vanishing minors of A satisfy the duality conditions of Eq. (14.51 ). Furthermore, any such A 
together with the phase parameters in the set K^, of Example 14.7 1 would generate a 3-soliton solution u{x,y,t) 
with soliton parameters p^. This solution is unique up to space-time translations corresponding to different 
choices for the phase constants 9^,0 1 m = 1 , . . . 6 in Eq. (11.51 ). 

We make a few observations from the above example. First, note that the rank conditions together 
with non-negativity of the maximal minors completely determine which maximal minors of A are zero, and 
which are non-zero. Secondly, the non-vanishing maximal minors are chosen to be positive by expressing 
them in terms of a suitable set of freely prescribed parameters which include certain matrix elements of 
A as well as certain combinations of the elements of A. This parametrization completely determines the 
coefficient matrix A which is in RREF and all the entries in its non-pivot columns can be expressed as 
appropriate maximal minors. Thirdly, both rank conditions (Propositions I2.6f i) and (ii)) are applied to the 
same index pair [e„,^„] since it labels a pair of asymptotic hne-solitons as \y\ 00. Instead, one could also 
apply the rank conditions from either Proposition 12. 6r i) or (ii) to each index pair, and recover the remaining 
information from the duality condition Eq. (14.51 ). For example, since columns 1 and 2 are proportional 
for the coefficient matrix A of Example 14. 8 [ the minors A(125) = A(126) = 0. Then the duality condition 
implies that A(346) = A(345) = as we found in Step 2 above. We also point out here that in Ref. [ HI, 
the A/^-soliton x-function was required to satisfy a set of conditions, which was referred to as the 'W-soliton 
condition" (Definition 4.2 of |14|). These conditions follow directly from the rank and duality conditions 
discussed in this article. We reiterate the above observations through another example. 

Example 4.9 In this example, we construct a 4-soliton solution starting with the following set of soliton 
parameters: 7^4 = = {(1, -3), (|, -|), (|, ^), (1,2)} G 5o/(4). 

Step 1. First, we construct the set 114 = {(~2, — 1), (— 1,0), (— ^, 1)(^, |)} from p4 using kf = ^{c„i: 
ttn). Then we obtain the ordered set = {— 2, — |, — 1, — ^,0, 5, 1, |} = {k\, . . . ^k%} by sorting the phase 
parameters in 114 which can now be re -expressed as 114 = {(^1,^3), (^21^5), (^41^7), (^61 ^s)}- This gives 
the correspondence 114 ^ (13)(25)(47)(68) = 71 G /s where, 7i is expressed as products of disjoint 2-cycles. 
The asymptotic line-solitons are identified by the set {[1,3], [2,5], [4,7], [6,8]} of index pairs. 

Step 2. We proceed to construct the coefficient matrix A that will generate the x-function of the 4-soliton 
solution. As in the previous example, we start with the 4x8 matrix A satisfying Condition |2.2[ 

-vi -w\\ 

-V3 -W3 ' 
V4 W4 J 

whose pivots and non-pivot indices are {^1,^2, ^3,^4} = {1,2,4,6}, and {^1,^2,^3,^4} = {3,5,7,8}, re- 
spectively, and where z\,Z2,u\,U2,uj and v,- , w, , / = 1 , . . . , 4 are non-negative reals to be determined. 



/I -zi Ml 

1 Z2 -M2 

1 M3 

VO 1 
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Figure 4.2: (a) a 3-soliton solution of Example 14. 81 at t = 12, generated by the coefficient matrix A with parameter 
values: z = V2 = W2 = 1, V3 = 3, vv'3 — ^ 

coefficient matrix A are given by zi =u\ = U2 = vi = V2 = vj, — I, Z2 — U3 ^ V4 ~ wi 



^; (b) a 4-soliton solution of Example 14.91 at t = 20. The parameters of the 



W2 = W3= 2, W4 = 3. 



Step 3. We then apply the rank conditions from Proposition 12.61 systematically to each soliton index pair 
and collect all the information regarding the unknown entries in terms of zero and non-zero minors of 
A. Here we indicate only the essential steps. 

Consider the index pair [1,3] and the associated sub-matrix X [13] which consists of columns 4, . . . , 8 of 
A above. Since rank(X [13]) < A^ - 1 = 3, the maximal minors: A(Zi, . .. ,/4) =0, {h,. . . ,14} C {4,5,6,7,8}. 
For the index pair [6, 8], the sub-matrix i'[68] = (— vi, V2, — V3, V4)^ is clearly of rank 1. Then from Proposi- 
tion l2.6f ii) rank(F[68]|6,8) = 2, which implies that A(Z,6,7,8) = for any / G [8]. Putting these conditions 
together, it can be shown that all vanishing minors of A obtained via Proposition 12.61 are generated by the 
following relations 



UiVi-UjVi = 0, 



UiW 



,.yj 



-UjWi = 0, ij= 1,2, ViWj-VjWi = 0, ij = 1,2,3. 



(4.6) 



Next we determine the non-zero entries of A. Note that zi 7^ since columns 2 and 3 of A are linearly 
independent. This is because the sub-matrix F[13] = (0, 1,0,0)^ associated with index pair [1,3], is of rank 
1, then Proposition I2.6f ii) implies that rank(F[13]|3) = 2. Consider now the sub-matrix F[47] whose rank 
is 2 since it is spanned by the linearly independent columns 5 and 6 of A. Then columns 4,5 and 6 are 
also linearly independent as rank(F[47]|4) = 3. Moreover, these columns also span the column space of the 
sub-matrix X[ 13] and form a basis which we denote by Pi 3. By similar reasoning, we find the basis sets P25 
and P47 for the column spaces of X[25] and X [47], respectively. These are given by 



)13 





U\ 







— U2 





1 


"3 





Vo 





1/ 



325 



/I 






-v,\ 

V2 

-"3 
V4 / 



347 



/I 






-Wl\ 

-W3 
W4 J 



Since both rank(X[13] 1 1) =rank(X[13]|3) =4, using the basis (3i3 we get A(1456) = M2 > 0, andA(3456) = 
Z2U\ — Z\U2 > where we also required that the non-zero maximal minors to be positive. In addition, if we 
take zi > (since zi / from above), then Z2 > 0, mi > 0, as well. Proceeding in a similar fashion with P25 
and P47 we find that that M3, V2, V3, V4, W3, W4 are all positive. Then from Eq. (14.61) . we finally conclude that 
all entries in the non-pivot columns of A are non-zero unlike the previous example. However, only 10 of 
the 13 non-zero elements of A are independent due to the constraints from Eq. (14.61 ). Furthermore, by direct 
computation using the matrix elements of A we find that all vanishing minors of A are in fact generated by 
the relations in Eq. (14.61 ) which arise from the rank conditions of Proposition 12.61 

As in Example |4. 81 it is possible to construct a set of positive parameters such that all non-zero maximal 
minors of A are polynomials in these parameters with positive coefficients. These are given by 



t\=z\., 



t2=Z2, 



Ui U2 
Zl Zl ' 



U2 

h = — ) 

Z2 



"3, 
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V2 V3 V3 W4 W3 W4 

t6 = ) = — , f8=V4, fg = , flO = 

U2 M3 M3 V4 M3 V4 



Figure |4.2f a) shows a 3-soliton solution with soliton parameters given in Example 14.71 and a coefficient 
matrix A constructed in Example 14.81 Figure I4.2f b) shows a 4-soliton solution with soliton parameters and 
matrix A from Example 14.91 It should be noted that the above examples are only illustrations of the general 
result asserting that for every pairing map % G I2N there exists a parametrized family of N x 2N matrices 
A which satisfy Condition 12.21 and generate an equivalence class of A/^-soliton solutions. The proof of this 
remarkable result will be given in a future work [61. 



4.4 Combinatorics of A^-soliton solutions 



Further refinement of the A/^-soliton classification scheme can be achieved by studying the combinatorial 
properties of the associated N x 2N coefficient matrix A. We have already shown that the A'^-soliton solution 
space is characterized by the set I2N of fixed point free involutions of the permutation group 52^- In turn, 
these involutions of S2N can be enumerated in terms of the various possible arrangements of the pivot and 
non-pivot columns of the A'^-soliton coefficient matrix A. A more geometric classification of the A/^-soliton so- 
lutions using Schubert decomposition of the real Grassmannian Gr{N, 2N) has been carried out in Ref. [ 14] 
where the arrangements of pivot and non-pivot indices were described in terms of Young diagrams. We 
resort to a more elementary treatment here and present the main results below. 

Proposition 4.10 Suppose that the index sets {ei, . • • ,^^1 and {g\,... ,g/v} with en < gn, n = \ , . ..N,form 
a disjoint partition of [2N], and label respectively, the pivot and non-pivot columns of a coefficient matrix A 
associated with a N-soliton solution ofKPII. Then, the following results hold. 

(i) The index set {ei,. . . jC^} is ordered as follows: I = ei < e2 < ■ ■ ■ < e^ < 2N. Moreover, the elements 
satisfy n <en <2n — I , n = \, . . . ,N. 

( ii) The total number of choices Cn for the ordered set {e\ , . . . , ca?} in item ( i) is given by the N-th Catalan 
number (see e.g., I[26\l }, 

( Hi) The set {gi ,. .. ,gN} is unordered, and the element g„ can be chosen in 2n — e„ ways for n= \ ,. .. ,N. 
Therefore, the number of possible choices for the unordered set {gi ,. .. ,gN} for each set {ei , e^} is 
given by 

N 

m(ei,...,e/v) = ]^(2?i-e„). (4.8) 



It follows from Proposition 14. 1 01 that the total number of distinct equivalence classes of A^-soliton solutions 
satisfies the curious combinatorial identity 

Fn:= £ m(^'i,...,e;v) = (2A^-l)!!, (4.9) 

n<e„<2n— 1 

which can be proved by following a similar line of argument that is provided below in the proof of Proposi- 
tion |4?T0tii). Items (i) and (iii) of Proposition 14. lOl were already proved in Ref. [14|. 

Proof. (Proposition \4.10U i)) Let En denote the set of all A/^-tuples {e\,...,eN) for which Proposition 
I4.10f i) holds. Then it is clear that \En\ = C^. Since ei = 1, each A/^-tuple contains one or more indices 
satisfying = 2m — I, m = I, . . . ,N. Then by sorting the elements of En according to the largest positive 
integer n € [A'^] such that Cn = 2n — I, we obtain the disjoint partition 

N 

En = \_\Wn where W„ = {(^i, . . . ,e/v) = 2?i — 1, Cm <2m— I, m> n} . 

n=\ 
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Note that W„ can be expressed as the direct product: W„ = En-i x {e„ = 2n — \} x E^^n, where the set 
£■„_! = {(^1 , . . . with j < ej < 27 — 1, 7 G [« — 1], is defined similarly as En, and where E!^^„ = 

{{e„+i,...,eN)} with 2?i + j- 1 < < 2(?i + 7) — 1, j£ [N — n]. If we define new indices e j :— en+j — 
{In — I), j = I, . . . ,N — n and relabel the elements of E^^,,, then it should be clear that 



En-u ^ {{ei,...,ej)\ j< ej < 2j - l} 



-N-n ■ 



Now the cardinalities of £"„-! and E!^^„ are C„_i and Ca?_„, respectively. Then it follows from above that 

|W;,| =C„-iCN-n, and 

N 

\En\ = Cn = Cn-lCN~n (4.10) 



n=l 



Eq. (14.101 ) gives a recursion relation for Cn, N > I with Co := 1. If C(z) = ^ Cnz^ is the generating 
function of the Cn, then Eq. (14.101 ) implies that C(z) satisfies zC^(z) — C(z) + 1 = 0, yielding 

^ 1-7T^ ^ - {2N)\ J, 

by choosing the root consistent with C(0) = Co = 1, then expanding it in power series. Finally, we obtain 
the desired result by equating the coefficients of the two power series for C(z). □ 

In view of Proposition 14. lOf i). it is natural to associate with each ordered set {ei , e^} a weight vector 
w and its length HwH defined by 

N 

w = {w\,W2,.--,wn), where w,, := e,j — « > 0, ne[N], and ||w|| = ^ (e„ — «) , (4.11) 

respectively. Note that the weights form a non-decreasing sequence: = w\ <W2<...<w/v and <Wn< 
n — I. Similarly, we associate the unordered set of non-pivot indices {g\ ,gf^} with an inversion vector a 
defined by 

a = (ai,a2---M. where On = \{gj\gj > gn, j <n}\, n£[N]. (4.12) 

The inversions satisfy < a„ < 2n — e^ — I, n = l,...,N. The upper limit of a„ follows from Proposi- 
tion K.lOf iii). by placing gn to the leftmost of the 2n — e„ available positions, and filling the remaining 
positions with indices gj such that j < n. Notice that the pair of vectors (w,a) are identical to the pair 
(y+,F^) of Young diagrams introduced in Ref |[T4l . 

The results of Proposition 14.101 together with the weight and inversion vectors provide a refinement of 
the classification scheme for the A/^-soliton solutions. We illustrate here the refined scheme for = 3. In this 
case, A is a 3 X 6 matrix with 3 pivots satisfying ei = 1, 2 < ^2 < 3, 3 < ^3 < 5. From Proposition 14. lOl ii). 
the total number of pivot configurations {^1,^21^3} is given by the Catalan number C3 = 5. Thus, there are 
5 subclasses of 3-soliton solutions depending on 5 distinct pivot configurations which are determined by 



w 


w 


{ei,e2,e3} 





(0,0,0) 


{1,2,3} 


1 


(0,0,1) 


{1,2,4} 


2 


(0,0,2) 


{1,2,5} 




(0,1,1) 


{1,3,4} 


3 


(0,1,2) 


{1,3,5} 



Table 1 : Weight vectors and possible pivot arrangements for 3-soliton solutions 
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{€1,62,63} 


{^1,^2,53} 


a 


3-soliton solution 


{1,2,3} 


{4,5,6} 
{4,6,5} 
{5,4,6} 
{5,6,4} 
{6,4,5} 
{6,5,4} 


(0,0,0) 
(0,0,1) 
(0,1,0) 
(0,0,2) 
(0,1,1) 
(0,1,2) 


{[1,4], [2,5], [3,6]} 
{[1,4], [2,6], [3,5]} 
{[1,5], [2,4], [3,6]} 
{[1,5], [2,6], [3,4]} 
{[1,6], [2,4], [3,5]} 
{[1,6], [2,5], [3,4]} 


{1,2,4} 


{3,5,6} 
{3,6,5} 
{5,3,6} 
{6,3,5} 


(0,0,0) 
(0,0,1) 
(0,1,0) 
(0,1,1) 


{[1,3], [2,5], [4,6]} 
{[1,3], [2,6], [4,5]} 
{[1,5], [2,3], [4,6]} 
{[1,6], [2,3], [4,5]} 


{1,2,5} 


{3,4,6} 
{4,3,6} 


(0,0,0) 
(0,1,0) 


{[1,3], [2,4], [5,6]} 
{[1,4], [2,3], [5,6]} 


{1,3,4} 


{2,5,6} 
{2,6,5} 


(0,0,0) 
(0,0,1) 


{[1,2], [3,5], [4,6]} 
{[1,2], [3,6], [4,5]} 


{1,3,5} 


{2,4,6} 


(0,0,0) 


{[1,2], [3,4], [5,6]} 



Table 2: The 15 distinct 3-soliton solutions 

the associated weight vectors. These subclasses are listed in Table[T]in increasing order of the length \\w\\ 
of the weight vector. For each of these pivot arrangements, Proposition I4.10f iii) gives the total number of 
distinct non-pivot configurations. Each set {^1,^2, ^3} in a given subclass, is distinguished by its unique 
inversion vector. The various 3-soliton solutions in each subclass is shown in Table |2]and the associated 
chord diagrams are presented in Figure 14.31 Note that instead of open chord diagrams, we use circles in 
Figure 1431 Here each straight chord replaces the pair of upper and lower chords connecting the same index 
pair [en,gn] in the self-dual, open chord diagrams of the A/^-soliton solutions (cf. Property [4.4f iv)). 

The subclass associated with the pivot configuration {1,2,3} is isomorphic to the permutation group 53 
acting on {4,5,6} to form the non-pivot sets {gi,g2,g3}- These are arranged in the second column of Table 
[2] according to the non-decreasing order of the L^-norm |a| = ai +02 + <73 of the respective inversion vectors 
in the second column of Table [2l The corresponding set of diagrams forming a hexagon in Figure 1431 is the 
permutahedron for 53. On the other hand when {61,62,63} = {1,2,4} the non-pivot index ^3 is chosen in 
2x3 — 4 = 2 ways; g2 is chosen in 2x2 — 2 = 2 ways; and obviously there is only one way to choose gi. 
Thus, there are only 4 (instead of 6) possible ways the non-pivot columns 3,5 and 6 can be arranged to form 
the set {g\,g2,g?,}- This is due to the restriction that ^3 7^ 3 because according to Proposition 14. 101 g3 must 
be greater than 63 = 4. The chord diagrams with the pivot set {1,2,4} form a square which is a subpolytope 
of the permutahedron of 53. We will discuss the poly tope structure for the A/^-soliton solutions in [6J. 

With the circular chord diagrams for the sets {6\ 6^} and {g\ ,...,gN}, one can easily find a ^-analog 
of the function m{6i, . . . ,6n) in Eq. (14.81) defined by 

N C,„ax 

m{6]_,...,6N){q) = W[2n-6n]q = Y^mcq'^ . 

Here rUc gives the number of the circular chord diagram having c crossings, and the maximum number of 
crossings for given {ei , . . . ,61^} is [6] 

N 

(-max — ^ • 
n=l 

For example, when {^1,^2, ^3} = {1,2,3} and {61,62,63} = {1,2,4}, we have 

m{\,2,3){q) =q^ + 2q^ + 2q+l and m{\,2,4){q) = q^ + 2q + I , 
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(436) 



(246) 




(346) 



(546) 



\ / 



(465) 



(356) 




(456) 

Figure 4.3: The closed chord diagrams for 3-soliton solutions. The dots on each diagram indicates the pivots 
{ei,e2;ej,), and the ordered numbers below the diagrams indicate (^1,^2,^3)- The number of the diagrams having 
the same number of crossings is given by the generating function Fj, {q) — + + 6^ + 5 where 5 is the Catalan 
number C3 — ^3(0). 



corresponding to the hexagon and square in Figure [431 Note that m{e\^. . . ,eN)ig = 1) = m{ei,. . . ^cm)- One 
can also define a ^-analog of the function Fj^ in Eq. ( 14.91 ) which gives the number of circular chord diagrams 
having c crossings for given N, i.e. the number of A/^-soliton solutions having c T-type interactions among N 
line solitons: 

Fn{(])-= m{ei,...,eN){q). 

{fi, 

For example, we have F3 (q) = + 3q^ + 6^ + 5 as is evident from counting the number of circular chord 
diagrams (from bottom to top) at each crossing level in Figure |431 Note also that Ffj{l) = Fj^ and Fiq{0) = 
Cn, the A^-th Catalan number which equals the number of possible pivot configurations [{ei as 
well as the number of chord diagrams with no crossings. The latter case yields the number of A/^-soliton 
solutions with only P- and 0-type interactions among the N line solitons. Furthermore, if we define the 
generating function by 

00 

F{q,x) := FN{q)x^ , 

N=0 

then it is possible to show that F{q,x) can be expressed by the continued fraction of the Stieltjes type. 



F{q,x) 



1 



1 



with Fo{q) = 1. A proof will be given in [6] (see also [|131 ). 
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